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A recent study revealed the dynamics of the charge sector of a one-dimensional quarter-filled elec- 
tronic system with extended Ifubbard interactions to be that of an effective pseudospin transverse- 
field Ising model (TFIM) in the strong coupling limit. With the twin motivations of studying 
the co-existing charge and spin order found in strongly correlated chain systems and the effects 
of inter-chain couplings, we investigate the phase diagram of coupled effective (TFIM) systems. 
A bosonisation and RG analysis for a two-leg TFIM ladder yields a rich phase diagram showing 
Wigner/Peierls charge order and Neel/dimer spin order. In a broad parameter regime, the orbital 
antiferromagnetic phase is found to be stable. An intermediate gapless phase of finite width is found 
to lie in between two charge-ordered gapped phases. Kosterlitz-Thouless transitions are found to 
lead from the gapless phase to either of the charge-ordered phases. A detailed analysis is also carried 
out for the dimensional crossover physics when many such pseudospin systems are coupled to one 
another. Importantly, the analysis reveals the key role of critical quantum fiuctuations in driving 
the strong dispersion in the transverse directions, as well as a T = deconfinement transition. Our 
work is potentially relevant for a unified description of a class of strongly correlated, quarter-filled 
chain and ladder systems. 



PACS numbers: 71.30.-|-h, 71.10.Pm 

I. INTRODUCTION 

Strongly correlated ladder systems are fascinating can- 
didates for studying the interplay of spin and charge or- 
dering, and their combined influence on the emergence 
of novel ground statesi. Several recent studies pro- 
vide experimental realisations of such systems, exhibit- 
ing diverse phenomena like charge order (CO), antiferro- 
magnetism (AF) and unconventional superconductivity 
(uSC) as functions of suitable control parameters-'- . The 
existence of several non-perturbative theoretical tech- 
niques in one dimension have also resulted in the study 
of the emergence of exotic phases from instabilities of the 
high-T Luttinger liquid^. Given that these systems are 
Mott insulators, longer-range Coulomb interactions are 
relevant in understanding CO/AF/uSC phases. Despite 
of some recent work^, a detailed understanding of the 
effects of longer-range interactions in quasi-lD systems 
is a largely unexplored problem. Further, attention has 
mostly focused on studies of models at 1/2-fiUing— . 

Here, we study an effective pseudospin model describ- 
ing charge degrees of freedom of a one-dimensional 1/4- 
filled electronic system. Such an effective model can 
be derived from an extended Hubbard model (i.e., with 
longer-range interactions) in a number of physically rel- 
evant cases. For example, it is reasonably well estab- 
lished that the structure of the material Na2V20^ is es- 
sentially that of a system of weakly coupled quarter-filled 
two-leg ladders — i^. While the spin sector is effectively 
quasi-one dimensional with each ladder corresponding to 
a spin chain, the charge sector is well described by a one- 
dimensional system of Ising pseudospins in a transverse 
field. The pseudospin degrees of freedom correspond to 
the position of a localised electron on a given rung of the 



two-leg ladder. 

Another relevant example is the Sri4^Cu24,04,i system, 
a ladder based material well described by a Hubbard-type 
fermionic model. Hitherto described by a Hubbard (or 
extended Hubbard) model at half-filling, recent experi- 
mental work by Abbamonte et al^ strongly suggests a 
very new scenario. Since pure Sri4Cu2404i is a Mott in- 
sulator with short-ranged antiferromagnetic correlations 
and a spin gap, nearest neighbour (nn) Coulomb inter- 
actions are in fact a necessary ingredient of a minimal 
model. Additionally, electronic charge-density-wave (e- 
CDW) order is inferred for zero doping, providing addi- 
tional evidence for inclusion of nn Coulomb interactions 
into a minimum effective model for these ladder systems. 
The e-CDW appears not to be driven by electron-phonon 
coupling, as shown by the absence of a detectable lattice 
distortion tied to the e-CDW. The relevance of inter- 
ladder coupling is also clearly revealed by this study. 

Similar arguments hold for the quasi-one dimensional 
charge transfer organics, which show e-CDW order in the 
Mott insulating phases in the generalised T — P phase 
diagram (substitution of different anions in the TMTSF- 
salts corresponds to varying chemical pressure)^*'. In 
studying a one-dimensional 1/4- filled fermionic model 
with extended interactions, it has been suggeste d^^i^^ 
that the same Ising pseudospin Hamiltonian may be a 
relevant starting point for the study of charge-order in or- 
ganics. Clearly, a half-filled Hubbard (or extended Hub- 
bard) model is inadequate when one seeks to understand 
Mott insulators with e-CDW (along with AF/dimerised) 
ground states. A quarter- filled, extended Hubbard model 
turns out to be a minimal model capable of describing 
such ground states — lii. 

While the weak coupling limit of the underlying ex- 
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tended Hubbard model has been recently studie d^^i^^ , we 
note that the real systems under consideration are gener- 
ically in the strong coupling regime of the model^i^i^. To 
the best of our knowledge, this regime has not been stud- 
ied in sufhcient detail. Another interest is to investigate 
the conditions under which short-coherence length su- 
perconductivity can arise by hole-doping a charge- and 
antiferromagnetically ordered Mott insulator. Again, 
this issue has been studied in sufficient detail only in 
the weak coupling limitii. Motivated by the above dis- 
cussion, we have recently studied a problem of coupled 
electronic chains^^, where each chain is described by 
an extended Hubbard model with a hopping term (of 
strength, t) and nearest- (Vi) and next-nearest neighbor 
(nnn) {V2) Coulomb interactions in addition to the lo- 
cal, Hubbard (U) interaction^?. Further, we studied the 
strong coupling regime, where U >> Vi,V2 >> t. In 
this regime, we derived an effective transverse-field Ising 
model (TFIM) in terms of pseudospins describing the 
charge degrees of freedom for a single chain — . In this 
work, our primary aim will be to study the effects of in- 
terchain couplings in this system of chains. As noted ear- 
lier, it has also been shown in Ref.(^) that a quarter-filled 
two-leg ladder in which the on-site Coulomb repulsion is 
the largest coupling can be mapped onto the same effec- 
tive TFIM Hamiltonian. Our results will, therefore, also 
be relevant to an understanding of such ladder systems. 

We begin in section II, by briefly reviewing for the 
sake of completeness, the derivation of the effective TFIM 
Hamiltonian for the charge sector of a system of 1 /4-filled 
electronic chains with strong extended Hubbard interac- 
tions and very weak nearest neighbour hopping. We refer 
readers to Ref.{— ) for a more detailed discussion. We 
then proceed with the effective pseudospin TFIM model 
for the charge degrees of freedom for a single chain^^ and 
couple neighboring chains to describe the physical situa- 
tions detailed above. In this way, a system of two effective 
pseudospin models coupled to one another is first anal- 
ysed in section II using abelian bosonisation and pertur- 
bative renormalisation group methods, to obtain a rich 
phase diagram showing different types of CO phases. At 
the same time, there remains open the intriguing possibil- 
ity that some of these gapped phases may themselves be 
separated from one another by non-trivial gapless phases 
of finite width^-. In what follows, we will encounter 
an example of this phenomenon in subsection IIB. The 
transitions from this gapless critical phase to either of 
the two charge-ordered phases are found to belong to 
the Kosterlitz-Thouless universality class. The question 
that naturally arises is whether such gapless, correlated, 
anistropic metallic phases survive when many such TFIM 
systems are coupled to one another so as to undergo a 
dimensional crossover. To answer this question reliably, 
we treat the inter-TFIM couplings at the level of the 
random phase approximation (RPA)^ '^^i^° in section HI. 
The finite T neighbourhood of the quantum critical point 
(QCP) is studied in detail. Interestingly, the inter-TFIM 
hopping coupling is seen to drive the system away from 



a gapped phase at T = to a gapless one at th QCP: 
critical quantum fluctuations drive the system through 
a deconfinement transition together with a dimensional 
crossover. In section IV, we present a comparison of our 
findings with some recent numerical works. Finally, we 
conclude in section V. 



II. TWO-LEG COUPLED TFIM LADDER 
MODEL 

We begin with the Hamiltonian for a spin chain system 

n 

+psis:,^2 + (1) 

where S':^, and S^^ and spin- 1/2 operators. The cou- 
plings {t, V, P) > are the nearest neighbour (nn) XY, 
the nearest neighbour Ising and the next nearest neigh- 
bour (nnn) Ising couplings respectively and h is the ex- 
ternal magnetic field. For h = 0, this Hamiltonian can 
be derived via a Jordan- Wigncr transformation on the 
charge degrees of freedom of a 1/ 4-filled extended Hub- 
bard model of electrons on a one-dimensional lattice in 
the strong-coupling limit: the on-site Hubbard interac- 
tion U ^ 00, while the nn (V) and nnn (P) density- 
density interaction couplings and the nn electron hopping 
(t) are all taken to be finite — 

H = ^[-^(c|ci+i -I- h.c) + VniUi+i + Pnini+2] ■ (2) 

n 

We study the problem in the limit of strong-coupling 
where V,P»t (but where (V - 2P) - 2t) 

Let us begin by studying the case oft = O^i (we will be 
studying eq.(IT]) for the case of /i = in all that follows). 
It is easy to see that for the case of ^ > 2P, the ground 
state of the system is given by a Neel-ordered antiferro- 
magnetic (AF) state with two degenerate ground-states 
given by 

\AFGS1) = 1... + - + - + - + -...) 
\AFGS2) = 1...- + - + - + - + ...) (3) 

where we signify S^^ — 1/2, —1/2 by + and — respectively 
and we have explicitly shown the spin configuration in the 
site nos. —3 < n < 4 in the ground states. In the original 
electronic Hamiltonian eq. ([2|) , this AF order corresponds 
to a Wigner charge-ordering (CO) in the ground-state. 
Similarly, for the case oiV < 2P, the ground state of the 
system is given by a dimer-ordered (2, 2) statei^ with four 
degenerate ground-states given by 

I22G51) = 1...- + + -- + + -...) 
I22G52) = 1...-- + + -- + + ...) 
I22G53) = 1... + -- + + --+...) 
|22G^4) = !... + + -- + + --...) 

(4) 
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where we signify 5*,^^ = 1/2, —1/2 by + and — respectively 
and we have exphcitly shown the spin configuration in 
the site nos. —2 < n < 5 in the ground states. In the 
original electronic Hamiltonian cq.([2]), this (2,2) order 
corresponds to a Pcicrls CO in the ground-state. 

The effect of the XY terms in the Hamiltonian ^ on 
these ground states is now considered. Let us start with 
noting the effect of a XY term on a single nn spin-pair 
on the 4 degenerate (2, 2) ground-states; for purposes of 
brevity, we will denote the entire t{S^S!ll_^_l + SnySf^^i) 
term simply as Thus, 



i°'i|22G51) 




+ -...) = -!...- + .. 


•> 


t°'i|22G52) 
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+ + ...)= 




t"'i|22G53) 




- + ...)-^|... + -.. 
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t°'i|22G51) 
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t-^'°\22GSl) 




.. + + ...)= 




t"i'°|22G52) 




..-+...> = -l...+-. 
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t-^^°\22GS'3) 




..--...)=0 




t"i-"|22G54) 




.. + -...) = -!...- + . 


..)(.5) 



In a similar manner, we study the action of the opera- 
tor ^"""+1 on the two degenerate ground states of the AF 
ordered configuration as 



f-^\AFGSl) 
t°'^\AFGS2) 
t-^-°\AFGSl) 
t-^'"\AFGS2) 



= r 



-l.Oi 



4l 



•>(6) 



(5f - 5f_i)/2, 



Defining bond-pseudospins 
S^S^_-^ and r^^ = S^_-^ (which can be rewritten in 
terms of bond-fermionic operators in the original elec- 
tronic Hamiltonian eq.([2) as = {ui — rii_i)/2, = 
c\ci-i and = Cic|_]^ respectively), we can write the 
four degenerate ground states of the (2, 2) ordered con- 
figuration in terms of these bond pseudospins as 



|22GS'l) 
|22GS'2) 
|22GS'3) 
|22GS'4) 



0-0 + 0.. 
+ - +.. 
+ 0-0.. 
- + - .. 



(7) 



where we have denoted r,^ = 1/2 as + and — —1/2 
as — and have explicitly shown the pseudospin configu- 
rations on the bond nos. < n < 4. We can clearly 
see from eq.(I7]) that these four ground-states break up 
into two pairs of doubly degenerate (AF) orderings of 
the pseudospins defined on the odd bonds (|22G51) and 



|22GS'3)) on the even bonds (|22G52) and |22GS'4)) re- 
spectively. It is also simple to sec from eq.© that the 
action of the operator t"^^'" (for the nearest neighbour 
pair of sites given by {n — l,n)) on these four ground 
states is to flip a pseudospin defined on the bond n (ly- 
ing in between the pair of sites (n — 1, n)) or to have no 
effect at all. 

We can now similarly see that the two-degenerate 
ground states of the AF ordered configuration can be 
written in terms of the bond-pseudospins defined above 
as 



\AFGS1) 
\AFGS2) 



(8) 



where we have explicitly shown the pseudospin config- 
urations on the bond nos. < n < 4. From eq.®, 
we see that the two degenerate ground-states have anti- 
ferromagnetic ordering of pseudospins on nn bonds; this 
can equally well be understood in terms of the ferromag- 
netic ordering of pseudospins on the odd bonds and on 
the even bonds separately. Further, from eq.®, we can 
see that the action of the operator (for the near- 

est neighbour pair of sites given by (n — 1, n)) on these 
two ground states is again to flip a pseudospin defined 
on the odd (even) bond n (lying in between the pair of 
sites (n — 1, n)) against a background of ferromagnetically 
ordered configuration of pseudospins defined on the odd 
(even) bonds. 

Thus, we can model these pseudospin ordered ground 
states dZ]),® as well as all possible pseudospin- flip ex- 
citations above them (as given by action of operators of 
the type i""^'" dS]),®) with the effective Hamihoniani^ 



H = - 



^ [2tr: + {V- 2P)<r^+2] 

n^odd 

[2t< + (V^-2P)«+2] 

n Get) en 

5][2t< + (T/-2P)«+2] , 



(9) 



where n is the bond index. 

This is just the Ising model in a transverse field, which 
is exactly solvabl o^^'^^ and has been studied extensively 
in lD-i2^. If {V — 2P) > 0, the ground state is ferromag- 
netically ordered in r^, i.e, it corresponds to a Wigner 
CDW. For {V — 2P) < 0, the Peierls dimer order results 
in the ground state. At {V — 2P) < 2t, the quantum dis- 
ordered phase has short-ranged pseudospin correlations, 
and is a charge "valence-bond" liquid. The quantum crit- 
ical point at {V — 2P) ~ 2t separating these phases is a 
deconfined phase with gapless pseudospin (r) excitations, 
and power-law fall-off in the pseudospin-pseudospin cor- 
relation functions. Correspondingly, the density-density 
correlation function has a power-law singular behavior 
at low energy, with an exponent a = 1/4 character- 
istic of the 2D Ising model at criticality. The gap in 
the pseudospin spectrum on either side of the critical 
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point is given by A,- = 2\V - 2P - 2t\. Further, the 
quantum-critical behaviour extends to temperatures as 
high as T ^ and undergoes finite-temperature 

crossovers to the two gapped phases at T ^ I^t|- The 
dynamics of the spin sector of a 1/ 4- filled electronic sys- 
tem in the limit of strong correlations was also studied 
in Ref.(-^); we restrict ourselves to a few brief remarks 
on the spin sector in this work and direct the reader to 
Ref.(J^) for more details. Henceforth, we will focus pri- 
marily on the effects of interchain couplings on the charge 
sector of such chain systems. 



A. Bosonisation and RG analysis 



each chain but commute between chains. It can indeed 
be checked that the bosonised expression for the various 
pseudospin operators satisfy the required spin commu- 
tation relation a^'^^ . In terms of these fermions, upon 
denoting the chains as a,b =t:ii we find an effective 
Hamiltonian for the ID Hubbard model with an equal- 
spin pairing term and an on-site spin-flip term 



H = E(^i«V',+i,a + h.c) - ^ Y.^^L'f'^-'' + ^•'^) 



Thus, we now proceed with the effective pseudospin 
Hamiltonian for the charge sector 



H 



chain 



[2t rf + {V- 2P) r^r^^^ 



(10) 



where the Ising pseudospin coupling V — 2P has been 
replaced by V for convenience in all that follows. Ro- 
tating the pseudospin axis ^ r^, — > — and in- 
troducing a bond fermion repulsion C^_l X)i a f)5^a("«,a ~ 
ni-\-i.a){ni_b ~ "-i+Lb) as well as a bond fermion transfer 
term U T.^,a,b^ai4,a^(^^+^■A■^clb,i(^^+l.b,i +^-<^) between 
two such chain systems described by the indices (a, b) (re- 
call the relations between the pseudospins and the bond 
fermions in the original electronic model given earlier), 
we have the effective Hamiltonian for the charge sector 
of the coupled system in terms of a pseudospin ladder 
model 

H = -^Y.'^2tTl, + V rl^T^^.J 

j,a,b^a 

where a, 6 = 1, 2 is the chain index. Having explored the 
strong coupling limit of U± >> y in an earlier work—, 
we will explore the weak coupling scenario of U± << V 
below. We note that in the spin sector, coupling between 
the two systems leads to a 5 = 1/2 Heisenberg ladder- 
type model, and has been studied extensively by several 
authors — li^. Interchain spin coupling turns out to be 
relevant, opening up a spin gap: in this case two possi- 
bilities are known to result. The ground state has either 
short-ranged antiferromagnetic correlations (RVB) with 
no magnetic long-range order (LRO), or it spontaneously 
breaks translation symmetry, leading to dimerisation. 

We now return to our analysis of the charge sector. 
To begin, we introduce fermionic operators 4'j,a on each 
chain via a Jordan- Wigner transformation of the pseu- 
dospins. We will then proceed to bosonise the theory. 
During this procedure, one has to take care that the 
spin commutation relations are maintained, i.e., the new 
fermionic operators have anticommutation relations on 



where t and t± are the in-chain and inter-chain hopping 
parameters respectively, /i = — 2t the chemical poten- 
tial of the effective model (and unrelated to the filling of 
the underlying electronic model), U = —U± the on-site 
(Hubbard) interaction coupling and V — V/A the pairing 
strength. Note that while we treat the parameters t and 
V as independent parameters for the sake of generality, 
t = V = V vti our original model (jlip . Bosonising in 
the usual way (in terms of the usual charge p =] + [ 
and spin a =t ~ i variables), we obtain the effective 
low-energy bosonic Hamiltonian 



27r 



dx[vpKpiTrUpix)y + -^{dx<Pp{x)) 



+ ^ J dxiv^K^inllAx))' + ^{dM^))'] 

U r U f 

+ — ^ / dx cos(V8(f>a) + -T-^ / dxcos(V8(j)p) 
2TTa J 27ra J 

\^ f \f2 f 

H / dxcos(V29p)cos(V20c) / dxdx(t>n(x) 

2-Ka J IT J 

+ — [ dxcos{^/2(j)^)cos{\/2e„) + — [ dxcoa{V8e^) 
ira J TTa J 

+ — [ dxcos(\/2(/)^)cos(%/26'p) (13) 
na J 



where H„ 



IdxOp 



VpKp 



v„K„ , Vp/Kp = vf{1 + and v^/K^ = ^^^(l - 

■:^^)- Among the various cosine potentials, we have 
the usual spin-flip backscattering cos(\/8</'cr) and Umk- 
lapp cos(-\/8(/>p) terms as well as the triplet supercon- 
ducting cos(-\/20p) cos(-\/2^cr) term. The chemical po- 
tential term can be absorbed by performing the shift 

4>p ^ 't>p + ^^""^ x. The cosine potentials with couplings 

Vi and V2 are generated under RG by the t± and V terms. 
Using the operator product expansion^, we find the RG 
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equations for the various couplings to second-order as 



B. Phase Diagram for repulsive inter-TFIM 
coupling 



(2 - 2K„)Ua -{^- K,)ti 



dl 

dy_ 

nr 

dtj_ 

IT 

dVi 

dl 



K, 



1 VV2 



(2- 



\2ti 



(14) 



where ah coupUngs have been normahsed with respect 
to the quantity and we have set 2'KVp — 1 for 

notational simphcity. The RG equations for the two in- 
teraction parameters [Kp, K„), the two velocities {vp, v„) 
as well as the parameter S — Kpfi/vp are found to be 



For repulsive interactions {U± > 0) between the bond- 
fermions, the a sector is massless and flows under RG 
to the fixed point value iT* > 1 and 1/2 < < 1. At 
1/2-filling (for the bond-fermions), the couplings Up, V, 
t±, Vi and V2 are all relevant while 11^ is irrelevant. The 
competition to reach strong-coupling first is, however, 
mainly between Up, t± and V. We show below the phase 
diagram as derived from this analysis. 



1.1 1 . 2 O O TTs 

A'p>(l/A'„,l/4(X„ + l/i^„)) 



0.9- 

o.s; 

^p_^ l/AV >K,> 1/4(1/ K„ + /K,) 



. 6 
. 5 



A'„ < (l/i{l/Kf, + 1/Ka), l/4{Ka + 1/Ka)) 



dK„ 

dl 
dKp 

dl 

dVg 

~dr 

dvp 
HT 
dS 
~dl 



-KliUl + Vi + 4) + +tl + V^ 
-KplU5{l)a) + + 
-VgKg{Ul + Vi + tl) + ^iVi +tl + V') 
-VpKpUp2i6{l)a) + ^{Vi + V^2^ 



(15) 



where 5(1) = 5e', a is a short-distance cut-off like the 
lattice spacing and Jq{x),Ji{x) are Bessel functions^. 
The various second order correction terms arise from 
{dx4'p/cr)'^ and {dxOp/a)"^ that are generated under the RG 
transformations^i2£i22. As discussed in Ref.(— ), the com- 
peting influences of the various couplings on the renor- 
malisations of the interaction parameters {Kp, K^) can 
cause their values to either grow or decrease. In turn, 
this affects drastically the scaling dimensions of the var- 
ious couplings and can lead to the system undergoing a 
transition from one type of ordered phase (in which a par- 
ticular coupling grows the fastest to strong coupling) to 
another, with the passage being through a gapless (crit- 
ical) phase. The existence of such a critical (gapless) 
region in coupling space, lying in between two ordered 
(massive) phases, is revealed in a subsequent analysis. 
Further, such a deconfined phase can be thought of as 
the quasi-lD analog of that analysed in detail later via a 
RPA treatment when dealing with many coupled chains. 



FIG. 1: The RG phase diagram in the {Kcr,Kp) plane for 
repulsive interchain interactions {U± < 0). The three regions 
Kp < (1/4(1/A-, + 1/Kp), 1/4{K„ + l/K^)), \/Ka > Kp > 
1/4(1/A'^ + 1/Kp) and Kp > {l/K„, l/4(Jf^ + l/K„)) give 
the values of (Ka,Kp) for which the couplings Up, t± and V 
respectively are the fastest to grow under RG. 



In the phase diagram in Fig.([T]), the three lines with 
intercepts at {Ko- = 1, Kp = 1), {K^r = 1, Kp = 0.64) and 
{Ka = 1, Kp = 1/2) are the relations Kcr — ^/Kp, Kp = 
l/4(l/i4:<, + 1/Kp) and Kp = l/i{K^ + 1/K„) respec- 
tively. The regions Kp < {1/A{1/K^ + 1/Kp), 1/A{K„ + 
1/K,)), 1/K, > Kp > 1/A{1/K, + 1/Kp) and Kp > 
{1/Ka, l/4:{Ka + l/K„ j) signify the values ofKp and K„ 
for which Up (in-chain Wigner charge-ordered Mott insu- 
lator) , t± and V (in-chain Peierls charge-ordered Mott in- 
sulator of preformed bond-fermion pairs) respectively are 
the fastest to reach strong-coupling. The RG equations 
for the coupling Up, the interaction parameter Kp and the 
incommensuration parameter S are familiar from the lit- 
erature on commensurate-incommensurate transitions^^. 
For temperatures T >> Kpfi, the finite chemical poten- 
tial (arising from the non-zero transverse-field strength 
t) is unable to quench the Umklapp scattering processes, 
allowing for the growth of Up to strong-coupling. For 
T << Kpfi, the finite chemical potential cuts off the RG 
flow of Up, freezing the Umklapp scattering processes. 

For the case of t± being the most relevant coupling, 
we find our RG equations to be a non-trivial gener- 
alisation of those derived in2^ for the model of two 
coupled spinless fermion chains without the intrachain 
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V pairing term. The resulting picture then describes 
strong interchain two-particle correlations between bond- 
fermions sharing a rung. Following the analysis outlined 
jj^j29^ we conclude that this phase is a novel insulat- 
ing phase characterised by a mass gap and delocalised 
bond-fcrmions on rungs, resembling the orbital antifer- 
romagnetic ground state found in the spinless two-chain 
problem. This matches our finding of an orbital anti- 
ferromagnetic ground state in the strongly-coupled lad- 
der with dominant antiferromagnetic rung-couplings in 
an earlier worki^. Away from 1/2- filling (for the bond- 
fermions), the competition is mainly between t±_ and V. 
For V reaching strong-coupling ahead of t±^ , the system 
is in a channel triplet-spin singlet superconducting phase 
with mobile intra-chain hole pairs; for tj_ reaching strong- 
coupling ahead of V, we are currently unable to describe 
in more detail the dominant instability away from the or- 
bital antifcrromagnetism like insulating phase away from 
quarter-filling. 



C. Gapless phase driven by inter-TFIM hopping 

Interestingly, while studies of ladder models have 
shown a plethora of charge and spin-ordered gapped 
phases — li^ii^, there remains the intriguing possibility 
that some of these gapped phases may themselves be sep- 
arated from one another by non-trivial gapless phases of 
finite width^*. In what follows, we provide an explicit 
realisation of this scenario. The RG equations p4p - (fT51) 
reveal the existence of a non-trivial fixed point (FP) for 
any value of {Kp, K^) lying the in ranges 1/2 < Kp < 1, 
1 < < 2 -I- -\/3 and which is perturbatively accessi- 
ble from the trivial weak-coupling FP. While perturba- 
tive RG is known to have its limitations in such cases, 
the finding of a non-trivial FP is nevertheless reliable as 
long as the values of the various couplings at the FP are 
small (as indicated above). We note that a similar non- 
trivial fixed point was found in a study of the anisotropic 
Heisenberg spin- 1/2 chain in a magnetic field^, where 
the authors derived a set of RG equations which were 
very similar to those found in2^. Here, the non-trivial 
FP is given by 

; — K* + 1 V,* 



2K* 



V* 



aK;{cK*-{K* + iyab) , V,* 




b V* 



where a = 2-{K; + 1/K;)/2 , b = 2- [l/K* + l/i^;)/2 
and c = 2 — {K* + \/K*)/2 . Further, we can safely 
make the approximation of the renormalisations of Kp 
and Kcr being small at this non-trivial FP The sys- 
tem is gapless at this non-trivial FP as well as at points 
which flow to it. The trivial FP has 6 unstable direc- 
tions {Up,V,tj_,Vi,V2 and S), 1 stable directions {Uc) 
and 2 marginal directions {Kp and K^-). The non-trivial 
FP has 5 unstable directions, 2 stable directions and two 



marginal directions. The presence of the two stable di- 
rections at the non-trivial FP indicates the existence of 
a two-dimensional surface of gapless theories in the five- 
dimensional (Ua, V, i±, Vi, V2) coupling space. This gap- 
less phase is the analog of the "Floating Phase" found in 
the phase diagram of the ID axial nnn Ising model2S. We 
present in Fig. ([2]) below a RG fiow phase diagram which 
is projected onto the (Vi, t±) plane (a similar RG flow di- 
agram is found for the case of the anisotropic Heisenberg 
model in a magnetic fleld h^-' in the (a, h) plane, where 
a is the anisotropy parameter). 




(0,0) 



FIG. 2; The RG phase diagram in the {Vi,t±) plane. The 
thick line characterises the set of points which flow to the 
intermediate fixed point {Vi ,t*j_) shown by the filled circle. 
The thin lines show all RG flows which flow towards strong- 
coupling in the two phases I and II, characterised by the 
relevant couplings V and t± respectively. 



The regions I and II characterise all RG flows which 
do not fiow to the intermediate fixed point at {Vi,t*^). 
The transitions from the non-trivial FP towards either 
of phases I and II fall under the universality class of 
the Kosterlitz-Thouless type^. In region I, Vi flows to 
strong-couphng while t± decays; for 1/2 < ifp < 1 and 
K„ > 1, we know from the above discussion that in this 
region, the coupling V will reach strong-coupling flrst. In 
region II, both t±_ and Vi grow under RG, with the cou- 
pling t± being the first to reach strong-coupling. Thus, 
the RG trajectory leading to the intermediate fixed point 
represents a gapless phase separating the two gapped, 
charge-ordered phases I and II characterised by the rel- 
evant couplings V and t± respectively. Nonperturbative 
insight on such critical phases is also gained in section 
IV, when we treat the case of many such TFIM systems 
coupled to one another using the RPA method. 
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D. Phase diagram for attractive inter-TFIM 
coupling 

For attractive interactions {U± < 0) between the bond- 
fermions, we can carry out a similar analysis. In this case, 
we can see that Kp > 1 while < 1. Then, from the 
RG equations given above, we can see that the Umklapp 
coupling Up and Vi are irrelevant while the couplings t±_ , 
C/ct, V and V2 are relevant. The competition to reach 
strong-coupling first is, however, mainly between Ucr, t± 
and V. We show below the phase diagram at 1/2-filling 
for the bond-fermions as derived from this analysis. 

1.1 1.2 1.3 1.4 1.5 

0-9^ \ > (l/Kp, 1/4(1//Cp + 1/K„)) 



0.8- 
0.7- 



0.6; 

0.5- 



1/Kp > K^> 1/VS 



< (l/4(l/ifp + 1/K^), 1/V3) 



FIG. 3: The RG phase diagram in the {Kp,Ka) plane for 
attractive interchain interactions {U± > 0). The three regions 
< {1/A(1/K„ + 1/Kp),l/Vi), l/Kp > > 1/73 and 
Ka > (1/Kp, 1/4.(1/ Kp + 1/K„)) give the values of {K„, Kp) 
for which the couplings Ua, t_L and V respectively are the 
fastest to grow under RG. 

In the phase diagram in Fig. ([3]), the three lines with 



intercepts at [Kp = l,Ka 



{Kp 



0.64) 



and [Kp = l,Ka = l/\/3) are the relations = ^/Kp, 
K„ = 1/4:{1/K„ + l/Kp) and K„ = 1/^3 respectively. 
The regions K„ < (l/4(l/if, + 1/^3), 1/Kp > 

> 1/V3 and > {l/Kp, l/A{l/ Kp + l/ K,)) signify 
the values of Kp and for which Ua (rung-dimer in- 
sulator with in-chain Wigner charge-ordering), and 
V (insulator with in-chain dimers and Peierls charge- 
ordering) respectively are the fastest to reach strong- 
coupling. This matches our finding of a ground state 
with in-chain Wigner charge order and rung-dimcrs in the 
strongly-coupled ladder with large ferromagnetic rung- 
couplings in an earlier work^^. Away from 1/2-filling 
(for the bond-fermions), depending on which of the three 
couplings t±, V and Ua is the first to reach strong- 
coupling, the system exists either as a superconductor 
with intra-chain hole pairs {V) or a superconductor with 
rung-singlet hole pairs [U^) or a phase reached by fol- 
lowing the dominant instability away from the orbital 
antiferromagnetism like insulating phase (tj_) but which 
we are currently unable to describe in greater detail. 



III. COUPLED TFIM SYSTEMS: QUANTUM 
CRITICALITY, DIMENSIONAL CROSSOVER 
AND DECONFINEMENT 

Having studied the rich phase diagram of a 2-leg cou- 
pled TFIM system in the preceeding section in consider- 
able detail, we now proceed to investigate the case when 
many such TFIM systems are coupled to one another. 
This is done with an effort towards gaining an under- 
standing of how such coupled systems undergo a dimen- 
sional crossover from nearly isolated quasi-one dimen- 
sional TIFM systems to a anisotropic strongly coupled 
system in higher dimensions. Dimensional crossover in 
coupled spin systems has been studied using renormali- 
sation group (RG) arguments''''^^ and the random phase 
approximation (RPA) for Ising chains^^i^^, Heisenberg 
chains^i^i^ and ladders in a magnetic field'^'* and lad- 
ders with frustration'^^. Having used perturbative RG 
arguments in exploring the phase diagram of the 2-leg 
system earlier, we now employ the RPA method to study 
the passage to higher dimensionality. This is in keep- 
ing with the fact that dimensional crossover is essentially 
a nonperturbative phenomenon ^. At the same time, 
while the mean-field-like approach of RPA is exact only 
in infinite dimensions (i.e., infinite coordination number), 
its application to the physics of coupled quasi-lD spin 
systems for small coordination numbers (i.e., lower di- 
mensions) has met with success ^S'Sr^ -[^ jg ^^^^ worth 
noting that while a naive mean-field treatment of single 
particle-hopping between fermionic chains is not possible 
as a single fermion operator has no well-defined classical 
limit we are able to treat two-particle hopping pro- 
cesses between our underlying chain systems^^ (or, lad- 
der systems^iii^) systems via RPA by working with an 
effective theory in terms of pseudospins (as evidenced by 
the inter chain bond fermion hopping processes studied 
earlier for the case of the 2-leg ladder). This is justi- 
fied because the presence of the large on-site Hubbard 
coupling in our underlying model makes the single parti- 
cle hopping irrelevant (in a RG sense) while two-particle 
processes (including hopping terms) can be crucial in de- 
termining the phase diagram^. A full treatment includ- 
ing single-particle hopping will require a treatment using 
chain-dynamical mean field theory (c-DMFT)'^* and will 
be the focus of a future work. In what follows, we will 
follow instead the RPA method outlined in RefJ^. 

We treat the dynamics of the coupled spin system for 
the two inter-TFIM couplings, J± and in eq. fTT]) 
given above using the RPA method in turn. Beginning 
with the coupling t±_ , this method involves computing the 
dynamical spin susceptibility x of coupled system in 
the disordered phase as 



X(w, k, k±) = [XiM^^ fc) - t±{k±)Y 



(17) 



in terms of the frequency uj, the longitudinal and trans- 
verse wave- vectors k and k± respectively. xiD is the 
dynamical spin susceptibility of a single TFIM system, 
to be calculated assuming incipient order along the 
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direction in pseudospin space 



(18) 

and the transverse coupling i_L(fcj_) ^ z±t^{ki_ = 0), for 
each TFIM system having a coordination number of z±_ . 
Then, a divergence in the dynamical pseudospin correla- 
tion function xi^^^^^i-) signifies an instability towards 
the formation of an ordered state. However, before set- 
ting out with the calculations, it is worth pausing to con- 
sider first the likely effects of a transverse coupling like 
t^. As discussed earlier, the phase diagram at T = 
and tj^ = is simple, with an ordered phase for 2t <V 
{g — {2t — V)/V < 0), a quantum disordered phase for 
2t > V {g > Q) and a quantum critical point at 2t = V 
((7 = 0). A finite t± will cause the ordered phase to be 
extended to finite temperatures (with a critical temper- 
ature Tc for the case of 2t = V), with a first order phase 
boundary ending at a new quantum critical point (QCP) 
gc = Ac/V at T = Q^K As for the simple TFIMi, there 
exists a "quantum critical" region just above the QCP 
and to the right of the ordered phase, with a crossover to 
the disordered phase at finite T. This is shown schemat- 
ically in the T — g phase diagram given below in Fig. ([¥]). 
The transition belongs to the 3D Ising universality class 
while the QCP to the 4D Ising universality class^. 




, (ilaplcss , ' 
(Juniitiini , 
('ritical , 
licgiou /■' 



Disordered 
Phase 



9c (QCP) 



FIG. 4: The T — g phase diagram for the case of many TFIM 
systems coupled by a transverse coupling (or Jx ) . The T = 
Ordered Phase of the uncoupled TFIM is now extended, 
with a phase boundary which has a value Tc for the model at 
g = Q and a T = quantum critical point (QCP) &t g = gc- 
The hashed region immediately to the right of the ordered 
phase and just above the QCP is the gapless quantum critical 
region (described in the text). The dashed line represents 
a finite T crossover from the quantum critical region to a 
Disordered phase. 

In the following, we compute, via the RPA method, the 
quantities gc dX T = 0, Tc for the case oi g = /S./V = 0, 



the shape of the phase boundary near the QCP, the dy- 
namical spin susceptibility x(<^j k±) at the QCP as well 
as the dispersion in the transverse directions for small t± 
and close to the QCP^^. We focus our attention mainly 
on, and in the neighbourhood of, the QCP in order to 
stress the role played by the critical quantum fluctuations 
in determining the physics of deconfinement and dimen- 
sional crossover in our system. In this, we are aided by 
the integrability and conformal invariance of the TFIM 
model for 2t = V, t± = oii2ii^; this allows us to exploit 
the nonperturbative results for a single TFIM system (as 
long as A = |2t— 1/| << 1), while deahng with the physics 
of the transverse couplings at a mean-field level. At the 
same time, the spectrum and dispersion deep inside the 
ordered phase also proves to be fascinating^^. Specifi- 
cally, it has been demonstrated that far from the tran- 
sition line, dispersion in the transverse directions is very 
weak (i.e., the spectrum nearly one-dimensional), and a 
hidden Eg symmetry of the underlying exactly solvable 
model^ gives rise to a spectrum of eight massive particles 
(three of which should be experimentally observable) . 

We begin by computing the critical value of the trans- 
verse coupling gc at T = 0. For this, we can use 
the expression for slightly off-critical xiD for the case 
when the mass of the single TFIM system is very small 
(to = Ac << 1). This, for small oj, is given by 



ZoV{Ac/V)^/^ 
Uj2 _ (vk)^ _ A2 



(19) 



where the velocity v = Va (where a is the lattice spac- 
ing) and Zq = 1.8437. Then, from the divergence of the 
susceptibility of the coupled system, x(w, fc, fc^) 



Xi^(w,fc) = zj_t^{k± = 0), 
we get, for the case oi uj — — k, 

Zq 

where we've dropped the argument of k± 
the sake of convenience. Thus, we get 



(20) 



(21) 



in for 



where the constant ci = Z, 



gc ~ ci(— ^) 

4/7 _ 



4/7 



(22) 



1.42. Precisely the same 



expression for the mass in the ordered phase and very 
close to the QCP, A = gcV, is also obtained by carry- 
ing out a self-consistent treatment of the effective mag- 
netic field, h ~ z±t±{T^) in the TFIM for a single chain. 
For this, one uses the slightly off-critical susceptibility 
for the ID TFIM given earlier (eqn. (fT9|l ) with the mass 
A replaced by A(l + [h/Vf) ii. From this result, the 
authors of Ref.(— ) concluded that the dispersion in the 
transverse directions in the ordered phase and close to 
the QCP is much stronger than that deep in the ordered 
phase. 
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To calculate the Tc for the case of g 
= 0, A: = 0) at finite T 



0, we use the 



Xioii^ 0, fc 



-7/4 



(23) 



where the constant C2 = sin(7r/8)B^(l/16, 7/8) and 
li{x,y) is the Euler Beta function. Thus, we find, by 
using eqn. 



T 

^ c 

v 



(24) 



where the constant C3 — ciji^ix) = 2.12. Further, we can 
also compute the susceptibility x for the coupled system 
at the QCP for small k± by using the relation for the 
slightly off-critical xiD given earlier, eqn. ([19]), together 



with the relation — gV^Vti_ in eqn. p?)) . giving 



x{uj,k,kA_) 



ZoVg. 



1/4 



uj^ — (vk)"^ — (v± ■ k±y 



(25) 



where \va^\^ = {ZoVgl^* /2)dH±ikj_ = 0)/dk']_ is gained 
by a Taylor expansion to second order—. Further, the 
shape of the phase boundary at low T can be determined 
by using the xiD of the TFIM at low 



Xioi^.k) = 



{u + i/r^f - {vkf 



(26) 



where = Yr^^^^ ^'^'^ dephasing time due to quan- 



tum fluctuations. 
egn. ipO)) gives 



Then, for cj = = fc in XiDj the 



A 



In T = In m + In A, 

T 



(27) 



where A = f (p^ - 1)^^^. The expression ^ given 
above has an approximate solutionis 

A 



ln(l/A) -lnln(l/A) 



(28) 



This relation gives us the shape of the phase boundary 
for low T and close to the QCP. In this way, we have de- 
rived the key features of the T — g phase diagram for 
the case of the t±_ transverse coupfing (Fig.(|3])) given 
above. We can alsocarry out an identical RPA calcu- 
lation for the other transverse coupling, Jx, for the case 
when the single chain system is critical, 2t — V,t^ — Q in 
egun. pij! ^^. This is because this theory is again known 
to be integrable^*' and falls in the same universality class 
as that of the TFIM. We can thus determine an identical 
set of relations for (i) the critical coupling at T = 0, 
(ii) the critical temperature for the case of g = 0, (iii) 
the susceptibility %, (iv) dispersion in the transverse di- 
rections as well as (v) the shape of the phase boundary 
close to the QCP to those obtained earlier, but with 
replaced by J± everywhere. In this way, we obtain es- 
sentially the same T — g phase diagram for the case of 



the J± transverse coupling, with the only difference be- 
ing the fact that the spectrum of the ordered phase is 
now obtained from the solution of the critical ID TFIM 
in a longitudinal fieldlA. 

We can see from the results given above that, as we 
move along the interior of the ordered phase towards the 
phase boundary, the excitation gaps decrease together 
with a gradual growth of the dispersion in the transverse 
directions. At the QCP, the spectrum is gapless and this 
is refiected in the fact that the dynamical susceptibility 
X at the QCP can have at most logarithmic corrections. 
At low T, directly in the region above the QCP, the spec- 
trum and dynamics are mainly governed by the quantum 
critical point while the thermal excitations are described 
by the associated continuum quantum field theory. By 
this, we mean that for energy scales given by a; >> fcsT, 
the system still resembles a quantum critical one in this 
region of the phase diagram while for lo « ksT, a relax- 
ation of the dynamics (given by t^^) caused by quantum 
fluctuations sets in. The rapid growth of the dispersion in 
the transverse directions close to the QCP is the physics 
of the dimensional crossover in this system, while the 
vanishing of all mass gaps at the QCP is the physics of 
deconfinement. These results echo our earlier flnding of 
the critical phase in the phase diagram of the 2-leg TFIM 
system by starting from the V ordered phase and increas- 
ing t± . While the integrability of the TFIM allows us to 
make considerable progress in computing various quanti- 
ties, the role of the QCP in the mechanism responsible for 
the dimensional crossover and deconfinement appear to 
be robust. These results lead us, therefore, to the conclu- 
sion that critical quantum fluctuations associated with a 
QCP can quite gcnerically cause such critical (gapless) 
phases to emerge in higher dimensions from the ordered 
(gapped) phases of lower dimensional systems when cou- 
pled to one another. 



IV. COMPARISON WITH RECENT 
NUMERICAL WORKS 

We present here a brief discussion of the relevance of 
our work to some numerical investigations that have been 
carried out on 2-leg ladder systems as well as coupled 
TFIM systems. Vojta et al4^ studied the problem of a 
strongly correlated electronic problem at 1/4-filling and 
with extended Hubbard interactions (keeping only a nn 
repulsion) using the DMRG method. The phase diagram 
they obtained contains several phases with charge and/or 
spin excitation gaps. While a comparison of our work 
with this study is hindered by the fact that the DMRG 
analysis docs not have the crucial element of the nun re- 
pulsion (V2 in our work). Fig. 2 of that work reveals that 
for the case oi U >> Vi > t, the authors indeed find 
a charge-ordered CDW state (i.e, the Wigner charge or- 
dered state of Ref.(~)) with an excitation gap in the spin 
sector as well. This is in conformity with our finding of a 
Wigner charge-ordered state with a spin gap for the case 



10 



of the Up coupling being the most relevant under RG. 
Further, the t±_ of that study corresponds to the single- 
particle hopping between the legs while our work has 
focused on the effects of two-particle hopping. Finally, 
with the on-site Hubbard coupling, [/, being the largest 
in the problem, we are unable to see any phase-separated 
state in our phase diagram (as observed by Vojta et al.). 

In a slightly earlier work, Riera et al.^^ studied the 
case of 1/4-filled chain/ladder Hubbard and t — J sys- 
tems, but which also include Holstein and/or Peierls-type 
couplings to the underlying lattice. Their findings reveal 
coexisting charge and spin orders in both chain and lad- 
der systems. Specifically, for the case of their chain sys- 
tem, by keeping only on-site and nn repulsion together 
with an on-site Holstein-type coupling of the electronic 
density to a phonon field, their phase diagram (Figs. 4) 
reveal separate phases with Wigner as well Peierls-type 
charge order. This is in keeping with our findings, but 
the origin of the Peierls order in the two cases are dif- 
ferent: in our work, it originates from the competition 
of the nnn coupling V2 with the nn coupling Vi , while in 
their work, it needs the Holstein coupling to the lattice. 
Riera et al. find a similar phase diagram (Figs. 5) for the 
case of an extended t — J model (i.e., including nnn t and 
J couplings) with a Holstein coupling. The addition of 
a Peierls-type coupling leads to a spin-Peierls instability, 
i.e., the formation of a dimerised spin order, which coex- 
ists with the Peierls-type charge order. Again, while this 
matches our findings, the origins are different. Qualita- 
tively similar conclusions are also reached by the authors 
in their study of an anisotropic 2-leg t — J ladder with an 
extended on-chain nn coupling and Holstein/Peierls-type 
lattice couphngs (Figs. 2 and 3). 

Finally, we comment on the very recent DMRG studies 
of Konik et al4^ on coupled TFIM systems in an effort 
at studying two dimensional coupled arrays of one di- 
mensional systems. By starting with a reliable spectrum 
truncation procedure for a single TFIM chain (which re- 
lies on the underlying continuum ID theory being either 
conformally invariant or gapped but integrable), the au- 
thors then implement an improvement of their DMRG 
algorithm using first-order perturbative RG arguments. 
Their results for a J± coupling of the TFIM chains con- 
firms the accuracy of the RPA analysis of Ref.(^^) and 
the present work in computing quantities like the single 
chain excitation gap (which is found to vanish at a critical 
Jj_) and dispersion of excitations in the coupled system 
as a function of J±. Their results indicate that the RPA 
method and the DMRG analysis agree very well upto val- 
ues of J± of the order of the gap. This method also ap- 
pears to give accurate values of critical exponents related 
to the ordering transition. Thus, this numerical approach 
appears to provide a confirmation of the interplay of the 
QCP in the TFIM and the transverse coupling in driving 
the dimensional crossover and deconfinement transition. 
Such an approach, therefore, holds much promise for the 
numerical investigations of such phenomena in systems 
with similar ingredients. 



V. CONCLUSIONS 

To conclude, we have studied a model of strongly cor- 
related coupled quasi-lD systems at 1/4-filling using an 
effective pseudospin TFIM model. Using a bosonisation 
and RG analysis for a 2-leg ladder model, we find two dif- 
ferent types of charge/spin ordered ground states at 1/4- 
filling. Transverse bond-fermion hopping is found to sta- 
bilise a new, gapped (insulating) phase characterised by 
interchain two-particle coherence of a type resembling or- 
bital antiferromagnetismiSiSS. The spin fluctuations are 
described by a 5 = 1/2 Heisenberg ladder- type model 
for all cases studied here: the spin excitations are always 
massive. Away from this filling, either intra- or inter- 
chain superconductivity in a gapped spin background is 
found to be the stable ground state. We also find the 
existence of an intermediate gapless phase lying in be- 
tween two gapped, charge-ordered phases (characterised 
by the relevant couplings V and t± respectively) in the 
RG phase diagram of our model. 

We followed this up with a RPA analysis of the case 
when many such effective TFIM systems are coupled. We 
find that a transverse bond-fermion (i.e., two particle) 
hopping coupling causes the ordered phase to extend to 
finite temperatures, with the phase boundary ending at 
a T = QCP. Interestingly, the critical quantum fluc- 
tuations at the QCP, together with the transverse cou- 
pling are found to drive a deconflnement transition 
(at T — 0) together with a dimensional crossover char- 
acterised by strong dispersion of the excitations of any 
single TFIM into the transverse directions. The gap- 
less higher dimensional system in the quantum critical 
region lying just above the QCP is in conformity with 
our flnding of a critical phase in the 2-leg ladder model 
driven by t± discussed above. Similar conclusions are 
also reached for the effects of the transverse coupling J± 
on effective TFIM chains which are critical. The mech- 
anism driving the dimensional crossover and deconflne- 
ment appear to be generic and lead us to believe that 
similar mechanism should exist for the case of other quan- 
tum critical systems in lower dimension when coupled to 
one another. Our present analyses are, however, espe- 
cially relevant to coupled chain systems like the TMTSF 
organic systems as well as coupled ladder systems like 
Sri4-3;Co-xCu2404i and a — NaV205 or — ^^00.33^205 
(a superconductor) which exhibits charge/spin long range 
order at a; = and superconductivity beyond under pres- 
sure and/or doping^i^. 
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